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Abstract 

We obtain the specific heat in the third constraint scenario for a canonical ensem- 
ble of a nonextensive extreme relativistic ideal gas in a closed form. The canonical 
ensemble of N particles in D dimensions is well-defined for the choice of the defor- 
mation parameter in the range < q < 1 + jjj^. For a relativistic ideal gas with 
particles of arbitrary mass a perturbative scheme in the nonextensivity parameter 
(1— q) is developed by employing an infinite product expansion of the g-exponential, 
and a direct transformation of the internal energy from the second to the third con- 
straint picture. All thermodynamic quantities may be uniformly evaluated to any 
desired perturbative order. 



PACS Number(s): 05.20.-y, 05.70.-a 
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I Introduction 



Tsallis pQ proposed nonextensive statistical mechanics by generalizing the functional form 
of the Boltzmann-Gibbs entropy as 

= k\n q W, qeR + , (1.1) 

where the deformation parameter q is taken to be a real positive number as this ensures [2] 
the stability of the Tsallis entropy. In (11. ip the quantity k is the generalized Boltzmann 
constant, and W denotes the weight. The entropy (II .ip satisfy a nonlinear, inhomogeneous 
relation 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B), (1.2) 

where A and B refer to statistically independent systems. The nonextensivity of the 
entropy manifest in (ll.2p is governed by the parameter (1 — q). The Boltzmann-Gibbs 
statistics is recovered in the q — > 1 limit. 

The nonextensive statistical mechanics has found wide-ranging applications in studies 
of the systems exhibiting long range interactions [3] , long time microscopic memory effects 
[I], anomalous diffusion |5j, nonequilibrium phenomena [6] and so on. For instance, the 
formation of a new hadronic state of matter known as the quark-gluon plasma that occurs 
in the early stage of the relativistic hadronic collisions exemplifies long range interactions 
as well as long time memory effects [7], and, consequently, the nonextensive statistical 
mechanics is expected to be more appropriate there than the classical Boltzmann-Gibbs 
statistics. The rapidity spectrum obtained by using the Tsallis distribution is found [8] 
to be in good agreement with the experimental data. The data on the distribution of 
transverse momentum of hadrons as well as the differential cross sections in high energy 
e + e~ collisions bear close resemblance with the theoretical analysis [HI CD] based on the 
Tsallis nonextensive statistical mechanics. In the context of many body systems endowed 
with self-gravitating long range interactions it has been observed [TT] that the power law 
distributions may be achieved using a ^-kinetic theory based on the Tsallis statistics. As a 
groundwork for complete theoretical understanding of the relativistic heavy ion collisions 
in the high energy physics regime, and also for possible applications in astrophysics, it 
is imperative to study the relativistic ideal gas in the context of nonextensive statistical 
mechanics. This investigation has been initiated in [12]. These authors observe that the 
grand canonical partition function exhibits an essential singularity for q > 1 region, and, 
consequently, they claim that, in the said region, the nonextensive relativistic ideal gas 
does not exist. Making a slight departure from the arguments in [12], we, in the present 
work, consider a canonical ensemble of a fixed number of N molecules of mass m of a 
relativistic ideal gas subject to the Tsallis statistics. We produce an exact evaluation of 
the canonical specific heat in the extreme relativistic case {m — > 0). In D dimensions 
the generalized partition function is nonsingular in the region < q < 1 + jj^. For an 
arbitrary mass m we obtain the generalized partition function [13], and the thermody- 
namic quantities in the second and the third constraint pictures as perturbative series in 
the nonextensivity parameter (1 — q). Towards this purpose, we, as a calculational tool, 
disentangle the (/-exponential (12.41) employing a technique developed in [H], and previ- 
ously used [15] to obtain thermodynamic quantities of a nonrelativistic ideal gas obeying 
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the Tsallis statistics. In addition, we employ a direct transformation linking the internal 
energies in the second and the third constraint pictures that allows us to evaluate the 
thermodynamic quantities uniformly to any arbitrary prescribed perturbative order in 
(1 — q). For the sake of simplicity we produce them till the order (1 — q) 2 . 

The plan of this article is as follows. The extreme or ultra relativistic gas is discussed in 
Sec. [Til This is followed by our perturbative evaluation of the thermodynamic quantities 
for the general case of a relativisitic ideal gas in Sec. IIHI We conclude in Sec. HV] 

II Ultra relativistic ideal gas 

The Hamiltonian of a relativistic ideal gas with particles possessing D-dimensional mo- 
menta pi (i = 1, • • • , N) reads 

JV 



i=l 



^) = E mc2 (\/ 1 + (£) » = 1*1. (2- 1 ) 



where c is the velocity of light. In the extreme relativistic case m « Pi the Hamiltonian 
( 12. ip reduces to 



N 



H(p) = cJ2p*- (2-2) 



i=i 



The generalized partition function in the third constraint [13J approach is given by 

zf{^N) = T7 ^ nm [ d DN x d DN p ex Pq (-P H{P \- U ® 1 , /?=-, (2.3) 



N\h DN J * ^ q \ r c( 3 ) J ' ' kT' 

where h refers to the elementary cell of the one-dimensional phase space, and the deformed 
g-exponential is the inverse of the q- logarithm introduced in (11.11) : 

ex Vq (z) = (l + (l-q)z)^. (2.4) 

The series expansion for the deformed exponential reads [16] 

(2.5) 



°° z n n 

exp (z) = V rT7 , [n] q = v — 

\ n W- l-(l-g)(n-l) 



We briefly remark here that the following g-derivative 

acts on the monomials x n to produce the g-number [n] g defined in (12.51) : 

D q (x)x n = [n] q x n ~\ (2.7) 
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Consequently, the g-exponentials represented by the infinite series (12. 5p are the eigenfunc- 
tions of the g-derivative introduced in (12. 6ft : 

D q (x) exp Jax) = aexp Jax). (2.8) 



Incidentally, other (/-derivatives were introduced in Ref.|17j. 

The ensemble probability for an arbitrary energy Ej in the third constraint framework 

leads [13] to the following sum of the q- weights, referred to in (12. 3ft : 

^^J2(pfW,V,N)Y=(zfy' 9 . (2.10) 

j 

For the extreme relativistic Hamiltonian (12. 2h the generalized partition function (12.31) 

reads 

zf\p,V,N) = g Zm=0 (P,V,N) (l + (l- q )JLup)^ +DN (c®)™ (2.11) 

where Z m= o(f3, V, N) is the classical partition function for the ultra relativistic gas in 
arbitrary dimension D 

Z m=0 ({3,V,N) = W{3- DN , (2.12) 
and the parameters W and Q are given by 

i (2V^v{D)\ N n r(^ + i) 



N\ \(ch) D T(%) ) ' * (i- q) DN T (_L- + DN + l y 

The generalized partition function (12.111) has simple poles for the values of the deformation 
parameter q — 1+ -, n = 1, ...,DN. The number of singularities equals the number of 
degrees of freedom of the system. The generalized partition function is, therefore, well 
defined in the interval < q < 1 + As the number of particles N increases, we observe 
that the poles accumulate towards q — 1, the limiting value where statistical mechanics 
becomes extensive. The internal energy is defined [13] via the escort probability as 

U®(P,V,N) = (cVy^ipfiP)) 9 Ej. (2.13) 

j 

Converting the above sum to the phase space integration a la (12.31) and employing (12.101 
12. lip we now get 

r(3) 

U®(J3,V,N) = DN—. (2.14) 

Applications of the relations (12. 10112.131) produce the explicit solution for the sum of the 
g-weights: 

1-9 l+(,l-q)DN (l-q)DN 

c ( 3 ) = (W Q) i-a- q )DN (1 + (1 - q)DN) p-i- { i- q)D N. (2.15) 
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Figure 1: Temperature dependence of specific heat for fixed q = 0.95 and various N 



Above equation in conjunction with ( 12.14ft now produces the internal energy as 



Uj} \(3, V, N) = DN (wgy-v-DDN (i + (I - q )DNy-v-«) DN 0-T=v=&m . (2.16) 



The specific heat defined as 



da 



(3) 



dT 



reads 



cf 

DNk 1 - (1 - q)DN 



^(3) 



(2.17) 



I l-q l+(l-q)DN (l-q)DN 

(WG) 1 -^ DN (1 + (1 - q)DN)i-(i-o° N p- i-(i-«)dn . (2.18) 



DNk. The be- 



The extensive limit of the specific heat is readily obtained: C q 
haviour of the specific heat with respect to the dimensionless scaled temperature 

(l-q)DN 

t 



yl/D \ i-(i- q )DN 

kT 



ch 



(2.19) 



for various values of N and q values are shown in the Figs. (CQ) and (j2J), respectively. 
Following [20] we now obtain the specific heat at the physically important thermodynamic 
limit N^oo,V^oo,p=y— » finite: 
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Figure 2: Temperature dependence of specific heat for fixed N{= 3) and various q 



which implies that the classical q — > 1 limit and the thermodynamic limit do not commute. 
As remarked in |2QJ the iV-independent negative specific heat realized in ( 12.201) for the 
extreme relativistic perfect gas may be of consequence in astrophysical problems. 

We have also calculated the specific heat in the third constraint picture using a detour 
[TB"] via the second constraint framework. The results obtained through the direct evalua- 
tion of the generalized partition function and the internal energy transformation method 
are identical. Further discussions will appear in Sec. IIHI 

III Relativistic ideal gas: molecules of arbitrary mass 

The relativistic ideal gas containing particles of arbitrary mass m, and described by the 
Hamiltonian ( 12.11) has possible applications in heavy ion collisions [TtfTO] in nuclear physics. 
Moreover, in the context of self-gravitating systems in astrophysics they are relevant [11] . 
In the opposite limits m » pi and m « Pi the general relativistic perfect gas reduces to 
the nonrelativistic case and the extreme relativistic case, respectively. For the relativistic 
ideal gas we obtain the thermodynamic quantities perturbatively by treating (1 — q) as 
the series parameter. In the present section we employ the second constraint approach 
[T3] as an intermediate step for evaluating the thermodynamic quantities. The physical 
variables obtained in the second constraint method are transformed [T3] to the respective 
quantities in the third constraint framework by introducing a fictitious temperature f3' 
that provides for the correspondence between the alternate constraints. The ensemble 



6 



probability and the partition function in the second constraint [IB] approach read 

P?W,V,N) = eXM J E3 \ Zf = ex P? (-^,). (3.1) 

The above probability is linked PS] to that of the third constraint approach as 

pf\0)=pf 090, (3.2) 
where the general transformation rule for the temperature reads 

l-(l-?)/3'^|i 

In the above equation the internal energy in the second constraint approach is defined 
PS] as 

^) = £(P?W^. (3-4) 

i 

Converting the sum over states in ( 13. ip to an integration over the phase space for the 
Hamiltonian (12.11) the partition function in the second constraint approach may be recast 
as 

Zf\(3, V, AO = | d™* cP^p exp 9 (-/3 ff(p)). (3.5) 

The integral over the phase-space in (13.51) may be performed exactly in the extreme 
relativistic case in the m — > limit: 

*«Cftv,AD . = ^H,(g,v,J0l ;i _,)ir^£ + gA0 ). < 3 - 6 > 



m=0 



where the classical partition function Z m=0 ((3, V, N) for the extreme relativistic perfect 
gas is given in (12.121) . In the general case with particles of arbitrary mass we proceed 
towards evaluating the integral (13.51) perturbatively. Assuming the partition function 
(13.51) to be well behaved in the neighborhood of q — 1, we follow [THfT5] to disentangle 
the g-exponential ( 12 Ah as an infinite product series of ordinary exponentials: 



exp fl (-/3 H(p)) = exp ( - ]T (1 g) " (TH n (p) ) = V{d p ) ex P (-/3 H(p)). (3.7) 

\ n=l 

In ( 13.71) the operator valued series T>(d@) reads 



where djg = f3 n -j^. The operator ( 13.81) links the partition function (13. 5ft in the second 
constraint approach with the classical Boltzmann-Gibbs partition function Z((3, V, N) as 

Z®(J3,V,N) = V{dp) Z(/3,V,N), Z(P,V,N) = ^W, V, 1)) N . (3.9) 
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Henceforth, for the purpose of simplicity, we consider D = 3. The phase-space integral of 
the classical partition function for a single particle reads 



Z(J3,V,1) = J d 3 xd 3 pexp (-£mc 2 |y 1 + {-^) ~ A \ > (3-10) 

and may be expressed [19] via the modified Bessel function of the second kind K n (z), n = 
2: 

t 7" \ 47rF exp(it) K 2 (u) h 9 . 

Z(P,V,1) = 7 , A = — , u = (3mc 2 . 3.11 

A d m mc 

For future use we express Z(f3, V, N) in a factorized form: 

ZCM^flO/W, M = ± ( >vM«>V> y. (3.12) 



A 3 y ' AH V « 

The recipe (13.91 13. lip now produces the partition function in the second constraint ap- 
proach as an infinite perturbative series in the nonextensivity parameter (1 — q): 

/ oo 2n \ 

Z^((3, V, AT) = Z(A V, iV) ( J>l) n (1 - q) n %(«) ) , (3-13) 

where we have defined 



v n=0 e=o 



K M = u irn- (3 - 14) 

The coefficients a n e(u) for the first few orders are listed below: 

a 00 (u) = 1, a w (u) = ^(1 + 3N)N -3N 2 u + ^(l + N)Nu 2 , 

a n (u) = -~(1-2N)N-N 2 u, a 12 (u) = ~(1 - N) N, 

a 20 (u) = 1{2 + 27N + 90N 2 + 81N 3 )N -^-{l + 3N)N 3 u 
8 2 

+ i(ll - 42N + 48N 2 + 54N 3 )Nu 2 + -(15 - 25N - 9N 2 )N 2 u 3 

— (2 - 3JV - 6N 2 - N 3 )Nu\ 
8 

a 2l (u) = --(11 -45N + 63N 2 -54N 3 )N- -(4- ION + 9N 2 )N 2 u 

+ _L(44 - 93N + 15N 2 + 54N 3 )A^ 2 + -(2 - 3N - N 2 )N 2 u 3 , 

a 22 (u) = --(25 - 63N + 56N 2 -18N 3 )N--(ll-20N + 9N 2 )N 2 u 
8 2 

+ ^(4-7N + 3N 3 )Nu 2 , 
a 23 (u) = -^(62 - 129N + 85N 2 -18N 3 ) N -^{2- 3N + N 2 )N 2 u, 
a 24 (u) = -1(6 - 11N + 6N 2 - N 3 ) N. (3.15) 



In the second constraint framework the internal energy (13 .4p may be recast as 



Uf\u) = -mc 2 -^- In, Zf\u). (3.16) 



The expressions (13.131) and (13.151) may now be employed to produce the following series 

(2) 

expansion for the internal energy Uq (u): 



Uf\u) f °° 



(Z(u, V, N)) 1 -* 5^(1 - q) n (jC(u)) , (3.17) 



v ra=0 



NkT 

where the first few coefficients read 

poo(u) = 3-«, p 01 (u) = l, p 10 {u) = -3Nu + ^(5-4N)u 2 + Nu 3 , 

Pll (u) = Q(2N-l)-5Nu + (l-N)u 2 , p 12 (u) = — ~ - - Nu, p 13 (u) = N - 1, 

p 20 (u) = -^(l + 3N)N + ^(l + 9N)Nu-^(ll-48N + 48N 2 )u 2 

— (26 - 41N)Nu 3 + -(28 - 48iV - 9N 2 )u 4 + (1 - N)Nu 5 , 
2 6 

P2i{u) = ^(22 - 102iV + 81iV 2 ) + ^(57- 102N)Nu - ^(163 - 3527V + 138iV 2 )w 2 

-3(6 - 7N)Nu 3 + (2 - 3N)u\ 

p 22 {u) = i(311 - 786iV + 516/V 2 ) + i(110-147/V)/V?i-i(175-312iV + 907V 2 )w 2 
4 2 6 

-4(1 - N)Nu 3 , 

p 23 (u) = ^( y 323 - 664N + 350N 2 ) + 3( y 8-9N)Nu-( y 5-8N + 2N 2 )u 2 , 

p 24 (u) = ^(ll-2(W + 9iV 2 ) + 3(l- AOA^, p 25 (u) = (3-5N + 2N 2 ). (3.18) 

In constructing the transformation leading to the thermodynamic quantities pertaining 
to the third constraint approach we first express the weight factor in (13.31) in terms of an 
integral over the phase space, and subsequently use a perturbative approach a la (13.91) : 

^' (Wv.^ ^^ 1 (3 ' 19) 

where the operator-valued series reads 

K{d p ) = 1 - (1 - q) (df + \df) + (1 - q) 2 [df + 5 -df + ±df) +.... (3.20) 

We follow (13.191) to compute the weight factor c^(/3), and subsequently substitute it in 
the transformation equation (13. 3ft . To eliminate a trivial kinematical dependence on the 
volume of the nonextensive system in our evaluation of its specific heat, we introduce an 
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appropriately scaled variable. Using the factorized form of the classical partition function 
(EED we define 

(3.21) 



u 



u 



and compute the inverse transformation in a series as 

/ oo 2n— 1 A 

/ = u 1 + £(1 - q) n Qni(u) (/C(u)) , (3.22) 
\ n=l e=o J 

where the first few perturbative coefficients read 

01o (u) = -(67V + ln/(u)) + 2iVu, g n (u) = 2N, 

2O (u) = ^(l + 97V + 21n/(u))iV+i(ln/(u)) 2 -3(5iV + ln/(u))iVu 

-^(9 + 5N - 4 In f(u))Nu 2 + 2N 2 u 3 , 

021 (u) = ^(7 + 22iV + 61n/(u))A^- 13A^ 2 u - 2(1 + A^)iVu 2 , 

022 (u) = ^(21 + 31iV + 41n/(u))iV-2iV 2 u, 023 (u) = 2(1 + N)N. (3.23) 

The internal energies in the second and the third constraint pictures may be directly 
interrelated. Employing the respective definitions (13 Ah and (12.131) in conjunction with the 
ensemble probabilities (13.11) and ( 12.91) pertaining to these two pictures, the internal energy 
in the third constraint scenario may be expressed in terms of the fictitious temperature 
P' as 

The above transformation method may be readily extended to any other thermodynamic 
average. The compendium of structures described in (I3.17H3.24"|) now produce the internal 
energy of an arbitrary relativistic gas in the third constraint picture as a perturbative series 
in (1 — q): 

n=0 \£=0 J 

To derive the above expression for the internal energy we have employed rescaling of the 
argument given in Eq. flA.ll) in the Appendix. The first few coefficients g n £ in (I3.25P are 



10 



listed below: 

£oo(» = 3-u, eoi(u) = l, 

eio(u) = 3N(3-u) + l(5 + 6A0u 2 -iVu 3 + (3 + u 2 )ln/(u), 

eil ( u ) = -6(l + A0 + 3iVu+(l + A0u 2 -31n/(u), 

M") = -y-6iV + iVu-ln/(u), £? 13 (u) = -(l + A0, 

g 20 {u) = -^(l-9A^)7V-18A^ 2 u-i(ll + 57iV-3iV 2 )u 2 

+1(17 + 12N)Nu 3 + 1(28 + 15JV + 6N 2 )u 4 - (1 + N)Nu 5 
2 6 

+ (3N(6 - u) - -(17 + 12iV)u 2 + 4iVu 3 + (1 + N)u 4 ) In f(u) 

+^(ln/(u)) 2 -2u 2 (ln/(u)) 2 , 

^ 21 (u) = 11 + 63N +^N 2 -6(3 + 2N)Nu 

-1(163 + 214iV + 114iV 2 )u 2 + 16(1 + AOiVu 3 + (2 + iV)u 4 

+ (6(3 + 2N) - (16 + 9iV)u - 16iVu 2 + 2iVu 3 ) In f(u) 

+1(9 -2u 2 ) (ln/(u)) 2 , 

g 22 (u) = 1(311 + 6007V + 2587V 2 ) - -(89 + 84/V)/Vu - -(175 + 162/V + 90iV 2 )u 2 
4 2 6 

+4(1 + N)Nu 3 + - (89 + 847V - 2(Wu - 8(1 + iV)u 2 ) In f(u) 
+5(ln/(u)) 2 , 

to(u) = -(323 + 4307V + 258/V 2 ) -22(1 + N)Nu- (5 + 4iV + 2/V 2 )u 2 

+(22(1 + N) - 2Nu) In /(u) + (In /(u)) 2 , 
£ 24 (u) = !(55 + 61iV + 40iV 2 ) -3(1 + N)Nu + 3{l + iV)ln/(u), 
Q25 (u) = 3 + 37V + 27V 2 . (3.26) 

With the internal energy (13. 25} I3.26P at hand, the definition of the specific heat (I2.17P in 
the third constraint scenario may be employed in conjunction with the scaling equation 
(I3.2ip for obtaining the following perturbative series in the nonextensivity parameter 

(!-<?): 



Cf(u) 



C (u) + (1 - q) d(u) + (1 - q) 2 C 2 (u) + . . . , (3.27) 



Nk 

where the extensive Boltzmann-Gibbs limit is given by Cq(u). We now enlist the first few 
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coefficients in the rhs of (13. 27ft : 
C (u) = 3 + u 2 -3/C(u) - (/C(u)) 2 , 

Ci(u) = 3iV(6-u)-^(17+12A^)u 2 + 4A^u 3 + (l + A^)u 4 + (3-4u 2 )ln/(u) 
3 ^6 + AN + 3 In / (u)) - 9Nu - 2 ^8 + 8N + In /(u)) u 2 + 2iVu 3 

[S( 



/C(u) 



+ -(89 + 84A^ + 201n/(u)) - lOTVu - 4(1 + iV)u ; 



(/C(u)) ; 



+2 ll(l + JV)+ln/(u; 



iVu 



(/C(u)) 3 + 3(l + iV)(/C(u))^ 



— (1 - 21iV)iV - 27iV 2 u + -(11 + 44iV - 8iV 2 )u 2 - ^(53 + 12iV)iVu 3 
--(219 + 232iV + 130iV 2 )u 4 + 19(1 + N)Nu 5 + (2 + A^)u 6 
+ ^ (5AN - 6Nu + (53 + 12iV)u 2 - 26iVu 3 - 38(1 + iV)u 4 + 4iVu 5 ) In /(u) 
+~ (3 + 13u 2 - 2u 4 ) (In /(u)) 2 - ^ (66 + 270iV - 81iV 2 + 12(9 + N) In f(u) 
+27(ln f(u)) 2 ) + 3(18 + 2N + 9 In f(u))Nu + - ( 1437 + 2044iV + 9067V 2 



+4(169 + 154iV) In /(u) + 60(ln /(u)) 2 Ju 2 - (169 + 154iV + 30 In f(u))Nu 3 

— (217 + 180iV + 93iV 2 + 24(1 + N) In /(u))u 4 + 8(1 + iV)iVu 5 /C(u) 
3 V / 

2221 + 1318 In /(u) + 158(ln /(u)) 2 + 2 f 1923 + 534 In /(u) + 636iv) N 



-2(659 + 534iV + 158 In /(u))JVu - 2 (l091 + 236 In /(u) + 8(ln f(u)f 
+4(1182 + 232 ln/(u) + 678^)^)^ + 4(118 + 8 ln/(u) + 116iV)iVu 3 



-4(17+ 137V + 6N 2 )u 4 



(/C(u)) 2 - ^ (4175 + 5488iV + 2922 TV 2 
+ (1324 + 1240iV) In /(u) + 84(ln /(u)) 2 ) - (331 + 3107V + 42 In f(u))Nu 
-- (700 + 672iV + 4027V 2 + 60(1 + N) In /(u))u 2 + 20(1 + N) Nu 3 (/C(u)) 3 



r 1 



- (2619 + 2996iV + 1850iV 2 + 4(111 + 109iV) In f(u) + 12 (In f(u)) 2 



-(111 + 109iV + 6 ln/(u))iVu - (30 + 27N + 16iV 2 )u 2 (/C(u)) 4 

167 + 1767V + 115iV 2 + 12(1 + N) In /(u) - 12(1 + N)Nu (/C(u)) E 
-5(3 + 3iV + 2iV 2 )(/C(u)) 6 . 



(3.28) 
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Figure 3: Dependence of specific heat on u for fixed q = 0.995 and various TV 



The dependence of the specific heat (13.281) on the dimensionless variable u for numerous 
values of N and q are displayed in the Figs. ([3]) and PJ, respectively. 

Above perturbative evaluation of the specific heat hinges on the series (I3.25P for the 
internal energy in the third constraint picture, and the thermodynamic transformation 
property ( 13. 24ft that is based on the equivalence of the ensemble probabilities ( 13.21) . An 
alternate determination of the specific heat utilizing the generalized partition function 
( 12. 3p is well-known [13]. Integration of the following thermodynamic relation 



duf ] d —(3) , . 

allows one to determine the internal energy in the third constraint picture, which, in 
turn, produces the specific heat via the definition (12.171) . In the present case, an explicit 
analytical integration of the differential equation ( 13.291) even as a perturbative series turns 
out to be difficult. The specific heat, however, may be directly extracted [15] from the 
generalized partition function as follows: 

^.-. ^f--' ,-,^ ,3,0) 
q ou 1 — q cm 1 — q 

even though an explicit evaluation of the internal energy Uf> by integrating (1X291) may 
not be feasible. The equivalence property in ( 13.301) is based on the equality ( 12.101) . A 
perturbative calculation of the specific heat employing (13.301) is now accomplished by using 
the sum of the (/-weights in the third constraint picture. The connection formula (13.21) 
for the probabilities along with the transformation equation (I3.22f) allows a systematic 
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Figure 4: Dependence of specific heat on u for fixed N = 3 changing q 



perturbative evaluation of the relevant quantity [Z ) x ~ q appearing in the rhs of (13.301) . A 
shortcoming of the method, however, is that it is imperative to evaluate the said quantity, 
say, at the order (1 — q) 3 for obtaining the specific heat, via (I3.30p . at the perturbative 
order (1 — q) 2 . 

Following the above description we now obtain a perturbative series for the generalized 
partition function raised to the exponent (1 — q) by using the equations (13.21 12.101 13.221 : 

(oo 2n-l \ 

^(l-?r^3n.(u)(/C(u))M. (3.31) 
n=0 1=0 / 
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The first few coefficients 3 n ^(u) in the above series read as follows: 
3oo(u) = 1, 3io(u) = JV(3 -u), 3u(u) = iV, 

32o(u) = ^(9iV-l)iV-6A^ 2 u+^(4 + 7iV)^u 2 -iV 2 u 3 + iV(3 + u 2 )ln/(u), 

32i (u) = -^(3 + 2A^)iV + 2^ 2 u+(l + ^)iVu 2 -3^1n/(u), 

3 22 (u) = ~(10 + 11N)N + N 2 u-N\nf(u), 3 23 (u) = -(1 + N)N, 

330 (u) = ^(l-27iV + 234A^ 2 )A^+|(l-21A^)A^ 2 u-i(33 + 168A^-72A^ 2 )A^u 2 
4 2 8 

+^(39 + 17N)N 2 u 3 + —(53 + 337V + 18iV 2 )iVu 4 - (1 + iV)iV 2 u 5 
6 12 

+ (3iV 2 (6 - u) - -(17 + 12A^)A^u 2 + 4iV 2 u 3 + (1 + AO^u 4 ) In /(u) 
-^u 2 (ln/(u)) 2 , 

3si (u) = i(33 + 174A^ + 18iV 2 )iV-^(3 + 2iV)A^ 2 u-^(445 + 600iV + 318iV 2 )iVu 2 

+ 15(1 + N)N 2 u 3 + (N + 2)iVu 4 + ^6(3 + 2N)N - 9N 2 u - 16(1 + iV)iVu 2 

+2iV 2 u 3 ) In /(u) + (6 - u 2 )iV(ln f(u)) 2 , 

332 (u) = ^(547+ 1050iV + 432jV 2 )iV- ^-(79 + 73iV)iV 2 u 
8 2 

— (169 + 159iV + 90iV 2 )Mi 2 + 4(1 + iVW 2 u 3 
6 

+ ( 89 + ^N)N - !0iV 2 u - 4(1 + N)Nu 4 ^ In f(u) + — N(ln /(u)) 2 , 

3 33 (u) = -^(907 + 1212iV+ 722iV 2 )iV- 21(1 + N)N 2 u- (5 + 47V + 2iV 2 )iVu 2 

+ (22(1 + N)N - 2iV 2 u) ln/(u) + iV(ln /(u)) 2 , 

3 34 (u) = i(107+119iV+78A^ 2 )A^-3(l + A^)A^ 2 u + 3(l + A^)iVuln/(u), 

3s 5 (u) = (3 + 3iV + 2iV 2 )iV. (3.32) 

The alternate method of evaluation of the specific heat may now be completed by sub- 
stituting the perturbative series (13.311 13.32p in the defining relation (I3.30p . The results 
obtained in this parallel procedure completely agrees with our previous evaluation pre- 
sented in (13.271) and ( 13.28ft . However, as we remarked earlier, an evaluation of the specific 
heat at the perturbative order (1 — q) n using the property ( I3.30f) necessitates evaluation of 
the relevant sum of the g-weights at the succeeding order (1 — q) n+l . Our prior evalu- 
ation of the specific heat based on direct determination (I3.25P of the internal energy that 
follows from the thermodynamic transformation property (13.241) requires computation of 
all pertinent quantities only up to the level (1 — q) n , and, therefore, has an advantage in 
this respect. 

We now obtain various physically relevant limits of the specific heat (I3.27p . For this 
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purpose we express the specific heat (13.271) in terms of the variable u introduced in (13.111) . 
rather than the scaled quantity u defined in ( 13.211) . The translation of our results (13.271 
I3.28f) to the corresponding quantities expressed in terms of the variable u may be most 
succinctly expressed by the functional form 

Cf\uMZ{u)) = Cf (u,ln/(u)), (3.33) 

where the classical partition function Z(u) is given in (13.121) . The lhs of (I3.33P is easily 
obtained by replacing the polynomials in u in the coefficients ( I3.28P by identical poly- 
nomials in the variable u, and simultaneous substitution of ln/(u) factors therein with 
the corresponding quantities in the variable In Z(u). For the sake of brevity we do not 
reproduce the full expression of the lhs in (13.331) . 

(i) For the nonrelativistic gas the limiting value is given by u ^> 1. We substitute the 
asymptotic expansion for the ratio of Bessel functions given in the Appendix (1A.2|) in the 



expression of the specific heat obtained via the replacement (13.331) in the equations (13.271 
I3.28H . We also note that in the m> 1 regime the classical relativistic partition function 
defined in ( 13.91) reduces pT9] to its nonrelativistic analog Z^ R (T, V, N): 



3N 

V N /?TrmkT\ ~ 

Z(T,V,N)U 1 = Z NR (T,V,N) = ^ (^=5-) . mi, 



where we have used the dimension D = 3. With the above facts in mind, we now obtain 
the nonrelativistic limit of the perturbative series for the specific heat of the ideal gas in 
the third constraint picture as 



2 u i 



NR = l + (l-q) (3N + \nZ NR ) 



3 Nk 

(i-q) 2 



(3N(2 - 21N) - 36N In Z NR - 4(ln Z NR ) 2 ^ + .... (3.35) 



The above series reproduces the specific heat of the nonrelativistic ideal gas [20] in the 
neighborhood of q — > 1. 

(ii) Another interesting limit is the extreme relativistic case that has been discussed 
in detail in Sec. [TT], where we have obtained the corresponding specific heat (12.181) in a 
closed form. Here, following the transition to the variable u that has been explained in 
the context of (13.331) . we consider the massless m — > limit of the perturbative series for 
the specific heat ( 13.271 13.28p . In the m — ► limit we observe that the classical partition 



function ( 13. 9ft reduces to 

Z(T,V,N)\ U ^ = Z m=0 (T,V,N), (3.36) 

where the classical partition function Z m= o(T, V, N) in the massless case is given in ( 12.121) . 
The u — > limiting value of the specific heat obtained from (13.271 13.281) reads 



l + (l-g)(6Ar + lnZ| m=0 ) 



37VA; 



(3N(1-21N) - 18iVlnZ| m=0 + (lnZ| m=0 ) 2 ) + . . . .(3.37) 
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As a further check on our results (13.271 13.28P the limiting series (13.371) in the vicinity of 
q — > 1 completely agrees with a perturbative expansion of the exact value ( 12. 18ft of the 
specific heat for the extreme relativistic ideal gas for the dimension D = 3. 

IV Remarks 

We have considered a canonical ensemble of N particles of a relativistic ideal gas, and 
found its specific heat in the third constraint scenario. In the extreme relativistic limit 
the generalized partition function, the internal energy, and, therefore, the specific heat 
may be exactly evaluated. This makes it possible for us to observe the nature of the 
singularities. The generalized partition function exhibits simple poles on the (/-plane at 
q = 1 H — , n = 1,2,..., DN, where the factor DN equals the number of degrees of 
freedom. The canonical ensemble is, consequently, well-defined in the parametric range 
< q < 1 + 5L. As N increases, the singularities approach the limit q — > 1 + that 
represents extensive statistical mechanics. This agrees with the earlier observation [12] 
that a grand canonical ensemble of a nonextensive relativistic ideal gas does not exist 
in the q > 1 regime. We also notice that in the present case of relativistic gas the 
thermodynamic limit and the extensive limit do not commute. For other systems similar 
results also hold [20lfT5] . 

For the general case of a relativistic ideal gas with arbitrary mass of the molecules we 
used a perturbative mechanism developed in [15J. The specific heat was obtained as a 
perturbative series in the nonextensivity parameter (1 — q) up to the second order. The 
evaluation was performed by using two different methods: a procedure that directly links 
the internal energies in the second and third constraint pictures via a transformation, 
and the traditional generalized partition function based approach. To our knowledge the 
former method has not been used earlier. The second approach requires the evaluation of 
the sum of the g-weights c*- 3 -* defined in ( I2.10p up to a perturbative order higher than the 
prescribed order of the specific heat, whereas in the first method based on the transfor- 
mation of internal energies all the quantities are uniformly computed up to the required 
perturbative order of evaluation of the specific heat. Both procedures generate identical 
perturbative series for the specific heat. As a bonus, the former approach easily produces 
the series for the internal energy in the third constraint picture. The known series of 
internal energy, in turn, produce the specific heat. The nonrelativistic and the extreme 
relativistic limits of the said perturbative series of the specific heat agree with known 
respective expressions. 

Using the integral representations of the gamma function, Prato [2T] observed that 
in the context of the second constraint picture, it is possible to connect the generalized 
partition function for a nonextensive statistical system with the partition function of the 
corresponding extensive (q — 1) system for both q > 1 and q < 1 domains. We now 
demonstrate such relations for the extreme relativistic ideal gas that has been exactly 
solved in our work. Moreover, such integral representations provide alternate derivation 
of the perturbative expansion scheme for the generalized partition function of a relativistic 
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ideal gas with particles of arbitrary mass. The Hilhorst integral for q > 1 region reads 

POO 

Z {2) (/3,V,N) = dt<5(t) Z(tp,V,N), (4.1) 
Jo 

where the kernel 0(£) is given by [21] 

0(f) = 1 — t^- 1 exp (--?-). (4.2) 

The classical partition function Z m= o(j3, V, N) for extreme relativistic ideal gas is given 
in (12.121) . The corresponding nonextensive partition function in the second constraint 
picture is obtained via (14.11) for the region q > 1: 



Z®(/3,V,N) = Z m=0 ((3,V,N) 



m=0 



r (^T ~ DN 

(q - i)™T(-y 



DN 1 

Z m=0 (/3, V, AT) | [ — -. (4.3) 

A J- 1 + (1 — 5) n 

71=1 V 



For an arbitrary mass of the molecules we consider the kernel (14.21) perturbatively. In the 
limit q — > 1 + , it behaves as a distribution comprising of delta function and its derivatives 
having support at f = 1: 

m = f(t -l) + I( 9 -1) (*-§*),(»-!) + .... (4.4) 

Substituting the above series of distributions in the Hilhorst integral (14. ip . we derive a 
perturbative connection formula between the generalized partition function in the second 
constraint scenario and the classical Boltzmann-Gibbs partition function. This connection 
formula agrees precisely with (13.91) obtained by using an infinite product expansion of the 
g-exponential. 

For the complementary q < 1 region, the generalized partition function in the second 
constraint framework is expressed [21] in terms of the Boltzmann-Gibbs partition function 
as a contour integral on a complex plane: 

Z® (P, V,N) = -^£dt 0(f) Z(tp, V, N), (4.5) 

where the kernel reads 

«(*) = ^l,^ exp ( - -*-) . (4.6) 

( g _!)-(— ) +1 V q-U 

The contour C on the complex plane in (14.5)) is comprised [21] of the segments {(oo, e), (f = 
e exp(zi?),0 < d < 2n), (e exp(z27r _ ), oo exp(i2n~))\ e — > 0}. Employing the kernel 
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(14. 6 p as before, we compute the generalized partition function for the extreme relativistic 
ideal gas in the q < 1 domain: 



ZW(I3,V,N) = Z m .„(J},V,N) 



m=0 



(l-q) DN T(jEi + DN S 



DN - 

As expected, the expressions (14.31) and (14.71) produce identical results for the general- 
ized partition function for the extreme relativistic ideal gas signalling continuity at the 
extensive parametric value q = 1. 



APPENDIX 

1. The translation of the argument of JC(u) defined in ( 13.14p is given by the following 
Taylor series: 

K{u) = /C(u) + (l-g)||(6iV + ln/(u))u 2 -2A^u 3 - (4(6iV + ln/(u)) - 8iVu 
-2Mi 2 ) /C(u) - ((14iV + In f{u)) - 2Nu) (/C(u)) 2 - 2iV(/C(u)) ; 

+ (1 - q) 2 ( - 3 Q(l + 69N + 22 In f(u))N + (In /(u)) 2 )u 2 + (757V 

+13 In /(u))iVu 3 + -(9 - 15N + AN In f{u))Nu 4 - 2iV 2 u 5 + (6(1 + 45iV 
+141n/(u))iV + 8(ln/(u)) 2 - 12(17iV + 3 In /(u))JVu - ^(57^ + 2307V 2 
+78N In f(u) + 2(ln /(u)) 2 )u 2 + (657V + 4 In f(u))Nu 3 - 2N(N - l)u 4 ) /C(u) 

+ (2 ( 87 + 975N + 258 ln f( u )) N + 6 ( ln /( u )) 2 - ( 247iV + 25 ln /( u )) ^ u 
-(23 + 38iV + 8 ln f{u))Nu 2 + 12iV 2 u 3 )/C(u) 2 + (-(105 + 574iV 
+ 102 ln/(u))iV + (ln /(u)) 2 - (89iV + 41n f(u))Nu - 4(1 + AOiVu 2 ) (/C(u)) 3 
+ (-(37 + 139iV + 12 ln f(u))N - 10iV 2 u) (/C(u)) 4 

+2(1 + 3iV)iV(/C(u)) 5 ) +.... (A.l) 
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2. The asymptotic expansion of the ratio of Bessel functions in the u >> 1 region reads 

Ki(u) ^ 3 15 15 135 45 7425 375 1095525 

K 2 (u) ~ ~ 2^ + 8^ ~ 8^ + i28n 4 + 32^ ~ 1024n 6 + 32m 7 ~~ 32768m 8 + " " ^ ' ^ 
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